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MAGNETOSTRICTION AND MAGNETOELASTIC QUANTUM OSCILLATIONS 

IN p-TYPE LEAD TELLU RIDE* 

byT. E. Thompson,! Paul R. Aron, B. S. Chandrasekhar, * and D. N. Langenberg** 

Lewis Research Center 

SUMMARY 

A detailed experimental and theoretical study of quantum oscillations in the mag- 
netostriction and Young's modulus of p-type lead telluride is presented. The valance 
band of lead telluride is approximated by a spheroidal, nonparabolic model in which the 
effects of strain on the valance band parameters are described by a deformation poten- 
tial model. Using appropriate thermodynamic derivatives of the modified Lifshitz- 
Kosevich expression for the oscillatory parts of the electronic free energy, it is shown 
that both types of oscillations arise mainly from relative shifts of the valance band max- 
ima due to shear strains, accompanied by intervalley charge transfer. Band param- 
eters derived from the periods, phases, and spin splitting of the oscillations are in 
generally good agreement with values reported by other workers. A detailed compari- 
son is made of the experimentally observed oscillation amplitudes with those predicted 
by theory, and satisfactory agreement is found. The ratio of the amplitudes of the two 
effects yields a value of the valance band deformation potential, E^ = 7. 9±1. 3 eV, in 
good agreement with a value found from piezoresistance experiments by Burke. 


INTRODUCTION 

The Landau quantization of the energy levels of the conduction electrons in a solid 
placed in a magnetic field forms the basis for some of the most powerful techniques for 
establishing the electronic s true time of solids (ref. 1) . One can construct, using the 
thermodynamic potential for the electronic system as calculated by Lifshitz and 
Kosevich (ref. 2), a free energy which is a function of the magnetic field, the tempera - 
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ture, and the elastic strain. The first derivatives of the free energy with respect to 
these three thermodynamic variables lead, respectively, to the magnetization, the en- 
tropy, and the elastic stress. Because of the quantization of the electronic states, the 
free energy is periodic in the reciprocal of the magnetic field. This oscillatory char- 
acter is consequently manifest in the magnetization, entropy, and stress. 

The oscillatory magnetization (de Haas -van Alphen effect) has been most extensively 
used in the study of Fermi surfaces in metals. More recently, the oscillatory entropy, 
which is revaled as temperature oscillations in a thermally isolated crystal subjected to 
a swept magnetic field, has also been developed as a useful tool for such studies. The 
third effect, which appears as an oscillatory magnetostriction (i.e. , the dimensional 
charge of a crystal placed in a magnetic field), forms one of the subjects of this paper. 
The elastic constants of the crystal, which are the second derivatives of the free energy 
with respect to strain, also show quantum oscillations, and form the second main topic 
of this paper . 

We have measured the oscillatory behavior of the magnetostriction and the Young’s 
modulus of single crystals of p-type lead telluride, a multivalley, degenerate semicon- 
ductor, and determined from these results the periods, phases, effective masses, and 
g-factors for this material. From the amplitudes of these two effects we have also de- 
termined the deformation potential for the valance band in lead telluride. The am- 
plitude analysis is considerably simplified by using the ratio of the amplitudes of the two 
effects. This method of obtaining deformation potentials in semiconductors, or the gen- 
eralized deformation parameters for Fermi surfaces in metals, is a new and important 
feature of the techniques described in this paper. 

Magnetostriction (MS) was first observed by Joule in 1842 in a bar of iron (ref. 3). 
The main body of work since that time has been on ferromagnetic materials in which the 
effect is largest: the strains in iron and nickel are typically 10”® to 10”® at the satura- 
tion field. The first observation of MS in a nonferromagnetic material was made by 
Kapitza in 1929 in a series of classic experiments in pulsed fields up to 250 kG (ref. 4). 
Kapitza observed a monotonic MS in bismuth, antimony, graphite, gallium, tin, and 
tungsten, finding the effect to be greatest in bismuth, where the strain e is ~10“® in a 
field of 15 tesla (150 kG) at room temperature. 

In 1963, Chandrasekhar pointed out that oscillatory MS should be observable in ma- 
terials that show the de Haas -van Alphen effect (ref. 5). Subsequently, Green and 
Chandrasekhar observed the first such oscillations in bismuth (ref. 6). Magnetostriction 
oscillations have since been observed in silver (ref. 7), arsenic (ref. 8), beryllium 
(ref. 9), bismuth (ref. 10), cadmium (ref. 11), chromium (ref. 12), copper (refs. 13 
and 6), gallium (P.A. Pentz private communication), gallium antimonide (ref. 14), lead 
telluride (refs. 15 and 16), antimony (ref. 17), lead (ref. 12 and private communications 
from M. M. Finkelstein and from A . J. Slavin), and zinc (ref. 18). One of us (P.R.A., 
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unpublished) has observed MS oscillations in gold. The oscillatory strain amplitudes 

— fi 

for these materials at 2. 5 tesla (25 kG) and temperatures near 4 K range between 10” 

Q -7 

and 10”. For lead telluride at this temperature and field, e ~ 10 . 

Magnetostriction in diamagnetic systems is related directly to the strain dependence 

of the Fermi surface. Shoenberg (see ref. 14) was one of the first to point out that the 

diamagnetic oscillatory magnetostriction could be analyzed in terms of the oscillatory 

magnetization M and the strain dependence of the appropriate extremal cross-sectional 

area of the Fermi surface: 
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where B is the magnetic field, s?. is a component of the elastic compliance tensor, 
and S m is the extremal cross-sectional area of the Fermi surface perpendicular to the 
applied magnetic field. (Throughout this report we use the standard "engineering" form 
of the reduced stress and strain tensors (ref. 19).) Equation (1) was derived using 
thermodynamic considerations together with the assumption that the strain dependence 
of the oscillatory diamagnetic free energy is dominated by the strain dependence of the 
S m which appears in the argument of the oscillatory factor in the free energy (see 
eq. (16)). The gallium antimonide and beryllium experiments were analyzed using equa- 
tion (1); a value for 3 In S m /3e. was found by measuring separately the magnetization 
and magnetostriction oscillations. In the present report we derive an expression for the 
magnetostriction in terms of 3 In S m /3e. directly from the total free energy, using a 
modified form of the Lifshitz-Kosevich expression for the diamagnetic free energy. We 
then use a simple multivalley, deformation potential model to relate 3 In S m /3 e. to the 
strain dependence of the lead telluride valance band. A similar method has been used 
by Aron and Chandrasekhar for bismuth (ref. 10) and by the Naval Ordnance Laboratory 
group for lead telluride (ref. 16). However, the results and conclusions in reference 16 
differ from ours at several points. An important new aspect of our analysis is the direct 
comparison of the MS oscillations with quantum oscillations in the elastic constant 
(ref. 20). We call these magnetoelastic (ME) oscillations to distinguish them from sound 
velocity oscillations which have also been observed (ref. 21). By taking the ratio of the 
amplitudes of these two effects, the analysis is considerably simplified in a manner sim- 
ilar to the simplification obtained by taking the ratio of the amplitudes of the MS oscilla- 
tions and. the magnetization oscillations in equation (1). 

Rodriguez (refs. 22 and 23) derived an expression for the velocity oscillations by 
considering the electronic contribution to the bulk modulus. He simply calculated the 
second derivative with respect to volume of the Lifshitz-Kosevich expression for the 
diamagnetic free energy of a free-electron gas, assuming the major contribution to this 
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derivative to come from the strain dependence of the S m , which appears in the oscil- 
latory factor, and using the free-electron expression 3 In S m /3V = -2/3. This approach 
is quite similar to those mentioned in the preceding section for calculating expressions 
for the MS oscillations. For the analysis of their beryllium experiments, Testardi and 
Condon (ref. 24) derived expressions relating elastic constant oscillations to magnetiza- 
tion oscillations. They measured separately the magnetization and velocity oscillations 
and from the combined data report values of 3 In S m /3ej. The ideas behind their argu- 
ments are quite similar to those used to derive equation (1) , where the magnetization 
and ME oscillations are related. Whereas Rodriguez's equations are limited to the free- 
electron gas, a ratio method such as that used by Testardi and Condon can be applied to 
an electronic system with a more general Fermi surface. Testardi and Condon also in- 
cluded B-H effects and the interaction between e and B in their analysis. 

In this report, we describe a spheroidal, nonparabolic, multivalley, deformation 
potential model of the valance band of lead telluride. We derive theoretical expressions 
in terms of our valance band model for the MS and ME oscillations from the appropriate 
strain derivatives of a modified form of the Lifshitz-Kosevich free energy. Band param- 
eters obtained from the experimental results are compared with results previously ob- 
tained by other workers . An analysis of the amplitudes of the MS and ME oscillations 
and their ratio is presented, leading to an experimental value of the valance band de- 
formation potential for PbTe. Some basic relations for the spheroidal, nonparabolic 
energy -band model are presented in the appendix. 


LEAD TELLURIDE VALANCE BAND 

Lead telluride (PbTe) is a cubic (NaCl structure , point group m3m) degenerate 

semiconductor characterized (at liquid-helium temperatures) by a small energy gap 

5 2 1 -1 

(~0.2 eV), high carrier mobilities 10 cm v sec ), and small effective masses 
(~0. 05 m Q ). A mass of experimental evidence (de Haas -van Alphen oscillations, 
Shubnikov-de Haas oscillations, Azbel-Kaner cyclotron resonance, interband magneto- 
optical absorption, and other optical and transport phenomena (see references in the be- 
ginning of ref. 25)) indicates that the Fermi surface of both p- and n-type PbTe is well 
described by (111) directed prolate spheroids centered on the L-points of the Brill ouin 
zone (fig. 1). There are eight half-spheroids or, equivalently, four full spheroids in 
the zone. 

Cuff, Ellett, and Kuglin (refs. 26 and 27) found a strong carrier concentration de- 
pendence of the (111) cyclotron mass for p-type PbTe, which indicates that the valance 
band of this material is nonparabolic . Dixon and Riedl (ref. 28) found similar results 
for the carrier concentration dependence of the electric -susceptibility hole mass of 
PbTe. In order to take nonparabolicity into account we use a model due largely to Kane 
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Figure 1. - The Brillouin zone of lead telluride showing 
the eight, (ill)- directed, half-spheroids that con- 
tain the carriers (either electrons or holes). 


(ref. 29) and Cohen (ref. 30). Considering the symmetry of PbTe and using a coordinate 
system where Pg is the momentum component along [111], and pj and pg are the 
two orthogonal components, Cohen’s equation (26) becomes 

2 2 / 2 
Pi + P2 _ ( E P 3 

2mj \ 2m^ 

where is an effective energy gap between the valance and conduction bands, m^ and 
mg are the band edge effective masses, and nig is the longitudinal effective mass of 
the conduction band, which is experimentally found to be of the same order of magnitude 
as mg (ref. 31). The Fermi surface is more easily visualized by writing equation (2) in 
the form 



E = 


2m t 


A. 

2mZ 


(3) 
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where 


and 


m t = m l 


+ ± + JL) 

E g 2E g m*y 


m l =m 3 


(4) 


(5) 


Equation (3) describes a surface of revolution about the 3 -axis, which has extremal 
cross sections at Pg = 0 and Pg = m^E g J(E/E g )(mg/m^ - 1) - lj . For the energy and 
mass values of our samples the second extremum does not exist, and the Fermi surface 
is spheroidal except for a slight bulging for large values of Pg. In this connection, we 
note that recent work by Schilz (ref. 32) indicates that when the carrier concentration be- 
comes large (p > 5x10 cm ), the valance band Fermi surface becomes cylindrical in 
shape. The extremal cross section perpendicular to the Pg axis, located at pg = 0, is 
given by 


S m3 = 2 ™1 E I 1 + 


( 6 ) 


The cyclotron orbits of importance in the present work do not traverse large values of 
Pg, so that the Pg dependence of the transverse mass indicated in equation (4) is small 
and can be neglected. Instead of equation (4) we therefore use 


m t = m i 



(7) 


Equations (3), (5), and (6) describe the spheroidal, nonparabolic model for the energy 
surfaces of p-type PbTe, which we use throughout the remainder of this report. As is 
conventional, we also define a mass anisotropy parameter K = m^/m^., which in our 

model has an energy dependence of the form K = (mg/m^l + (E/E g )J“*. 

A test of the validity of the model and an evaluation of the numerical values of the 
model parameters appropriate to p-PbTe can be made using our data and the data of 
Cuff, Ellett, and Kuglin (ref. 26). When the magnetic field is parallel to the unique axis 
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of the spheroid, the quantum oscillation period Pg is given by 


P 



eh 


cE F m i 



and the corresponding normalized cyclotron mass pg is 


Pg S 


1 

27rm o 



( 8 ) 


(9) 


(For derivations of these relations, see the appendix. In this report p is always used 
to denote cyclotron masses normalized to m Q , the free electron mass, and m is always 
used to denote unnormalized band masses. Notice that pgm Q / m^., ) Solving for Ep 
from equation (8), substituting this value of Ep into equation (9), and squaring the re- 
sult give the relation 




+ Cm. (P,E m)’ 1 
1' o g o 


( 10 ) 


-4 2 

where C = 4 eh/(m c) = 4.63x10“ electron volt per tesla. Thus, a plot of pf against 

_ U t 

Pg gives the band edge mass m^ and the energy gap E g . Such an analysis was car- 
ried out in reference 27 , and in figure 2 shows a plot of their data along with values 
found from this experiment. A least-squares fit of the data gives n^ = (0. 018±0. 003) 
m Q and E^ = 0.10±0.02 electron volt. Cuff, Ellett, Kuglin, and Williams (ref. 27) find 
a different value for m^. This discrepancy appears to originate from the fact that their 
equation (2) is not equivalent to our equation (10), although it appears that the same 
model is being considered. A similar study indicates that m, is essentially constant 

17 ^ 1 Q 

over the range of carrier concentrations 3.6X10 per cubic centimeter < p < 3. 5x10 
per cubic centimeter and has a value m^ = mg = (0.10±0.02) m Q . 

The direct gap A between the valance and conduction bands of PbTe at T = 4. 2 K 
has been determined in two different experiments. Mitchell, Palik, and Zemel (ref. 33) 
find A = 0. 190±0.002 electron volt from magneto-optical studies of epitaxial films and 
Butler and Calawa (ref. 34) find A = 0. 187 electron volt from magneto-emission studies 
with PbTe diode lasers. The fact that the effective gap E^ in our model is smaller than 
the measured direct optical gap A indicates that interactions between the valence band 
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Figure 2. - Square of normalized cyclotron mass plotted against 
reciprocal period for B|| dll]. The intercept and slope are 
used to calculated and Eg from equation (10). 

and the higher conduction band states are important in determining the valance band 
parameters. 

The interaction of the electron spins with the magnetic field splits each Landau level 
by an amount ±gehB/4 m Q c, where the electron g-factor depends on the electron spin- 
orbit interaction. The two series of quantum oscillations associated with the two spin 
orientations interfere, producing a single set of quantum oscillations with an amplitude 
controlled by the factor cos( 7 rgp /2) , where p is the normalized cyclotron mass appro- 
priate to the applied field direction. Knowledge of this factor is important in our ampli- 
tude analysis. 

The g-factor is related to the components of a second-rank tensor and in PbTe has a 
directional dependence of the form (refs. 35 and 36) 

g 2 = g 2 cos 2 # + g 2 sin 2 # (11) 

where g| ( and g^ are the g-factors for the field parallel and perpendicular to the [111] 
axis. The results of a number of different experiments (table I) indicate that at the band 
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TABLE I. - SPIN SPLITTING FACTORS AND g-VALUES FOR p-TYPE LEAD TELLURIDE 


Experiment 

Carrier 


Effective mass g-factor 


Band -edge 

g» 

source 

concentration. 





effective mass, 



cm -3 

(gp /2 )<too) 

^(100) 

g <100) 

g"P t /2 

** 3 


Ref. 34 

a n 



b 29 

















51*8 








This report; 

.9xl0 18 

0. 65±0. 05 

c 0. 04l±0. 004 

32*4 


0. 027*0. 002 

d 48±5 

sample X6 








Ref. 37 1 

3.0 

. 58±0. 01 

. 051*0.008 

23*5 | 

0.58*0.01 

. 036*0. 002 | 

32±2 

This report; 1 

3.5 

. 55*0.05 

c . 055±0. 002 

20±2 


. 038±0. 001 

d 29±3 

sample X5 








Ref. 32 

4 



— 

e . 58*0.03 

. 040*0.008 

29±6 

Ref. 32 
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— 

e . 53±0. 02 

.044*0.009 

24±5 


a p-n junction. 

b Assuming conduction and valence band g-factors equal. 
c Computed from eq. (A9). 

d Assuming gp to be the same along (111) and (100) . 
e Using a different splitting scheme than Schilz. 


edge, gjj s 50 and that g |( decreases with increasing carrier concentration. Burke, 
Houston, and Savage (ref. 37) find g^ = 7±2. The carrier concentration dependence of 
gn is related to the nonparabolicity of the bands; the interband interactions which deter- 
mine the g-factor are quite similar to those which determine the inverse effective mass. 
To a first approximation, both experimentally and theoretically (refs. 35 and 36) g|j 
varies with carrier concentration as 1/m^. Schilz (ref. 32) found the spin factor gp/2 
to be directionally independent. This is borne out by other work; for example, using 
the values of Burke, Houston, and Savage (ref. 37) and equations (A9) and (11) we 
compute gp/2 to be 0. 56 for the ellipsoid with the smallest extremal area when B is 
along [110]. This value is sensibly the same as the values 0. 58 shown in table I for the 
[100] and [111] directions. Since we were able to measure the spin splitting only with 
Bll[100], we will assume the amplitude factor cos(irgp/2) to be the same for the [100], 
[110], and [111] directions. 

The effect of strain on nondegenerate * energy levels in semiconductors can be de- 
scribed by a deformation potential tensor E (refs. 38 to 40). In this scheme the shift in 
the energy of the point v in a band is described to first order in the strain IT by 
5E V = E v : V! In general, each point in the band will have nine deformation potentials 


^There is, in fact, the spin degeneracy; however, for the rock salt structure, this 
degeneracy is never destroyed by strain deformation (ref. 41). 
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associated with it. However, this number is usually greatly reduced at a symmetry 
point. At the L-points in PbTe there are only two independent deformation potentials for 
each band, and we can express E v , using the notation of Herring (refs. 39 and 40), as 
E y = *d 5^ + E u u^uj where 6^ is the Kronecker delta and the unit vector u points 
to a particular L-point. This leads to an energy shift of the form 

H d + l S u) (e l + e 2 + e 3 } + “ H u (± e 4 ± e 5 ± e 6 } (12) 

where the signs depend upon u. The PbTe deformation potentials have been calculated 
by Ferreira (ref. 41). However, there has been only a limited amount of experimental 
work on determining the PbTe deformation potentials. Ferreira's results are shown in 
table II along with the deformation potentials derived from piezoresistance measurements 

TABLE II. - DEFORMATION POTENTIALS OF 
LEAD TELLURIDE 



Type of data 

Source 

Valence band 

Conduction band 



Deformation potential, eV 




~v 

~c 

~c 



~d 

~u 

~d 

~u 

Theory 

Ref. 42 

-8.93 

10.46 

-4.36 

8.29 

Experiment 

Ref. 43 


5.6 


3.8 




8.5 


— 


by Ilisavskii (ref. 42) and Burke (ref. 43). In order to find the deformation potential 
from piezoresistance, the elastic constants and the ratio of the mobilities perpendicular 
and parallel to the [111] axis must be known. Ilisavskii did not have good values for 
these two parameters, so we have recalculated his deformation potentials using Burke's 
latest values. The variation of the energy gap with pressure has been measured optically 
(ref. 44) giving ^3E^ + E^ - ^3E^ + = 9. 5 to 11. 4 electron volts, where the super- 

scripts refer to the valance or conduction bands. Ferreira's theoretical values give 
11.6 electron volts, in reasonable agreement with experiment. 

The band edge masses m^ and m^ are found in the work of Lin and Kleinman 
(ref. 25) to be approximately given by nig/m^ = M^/A, where the ]VL(i = 1,3) are 
appropriate momentum matrix elements. Since we expect these matrix elements to be 
insensitive to small strains, we find, using equation (12), 
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( 13 ) 


where 


3 In m 


1 , 3 _ 3 In A 


Se i 


3 e i 


G-i 

— for i = 1, 2, or 3 
A 




± — for i = 4, 5, or 6 
A 





(14) 


and the sign depends on the particular L -point considered. Since the g-values depend on 
the energy gap in a similar manner, they should have a similar type of strain dependence. 


THEORY OF OSCILLATING MECHANICAL EFFECTS 
Thermodynamic Relations 

We obtain theoretical expressions (e.g. , ref. 45) for the magnetostriction and mag- 
netoelastic oscillations by considering the free -energy density of the crystal in a mag- 
netic field as the sum of the zero-field elastic energy and the free energy of the diamag- 
netic electronic system £2: 


6 

F(e,T,H) = i > cf.e.e. + n(c,T,B) + 2wM 2 (15) 

2 Z-f 1] 1 ] 

M=1 

where the c?. are zero -field components of the elastic stiffness tensor. The magnetiza- 
tion can be obtained from this free energy as (ref. 45) M = -3F/3H = -3S2/3B and the 
stress as o\ = 3F/3e. The explicit dependence of the chemical potential on the strain, 
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temperature, and magnetic field is neglected since far from the quantum limit (the ap- 
proximation considered here) its contribution to the oscillatory free energy will be 
small. Since we are interested in oscillatory effects, we need consider only the oscilla- 
tory part of to. For this we use the Lifshitz-Kosevich expression in the modified form 
given by Williamson, Foner, and White (ref. 46): 



-*x 


(16) 


where A^ is a temperature -dependent amplitude function that varies slowly with field, 


- 2kT 


© 


3/2 


aV 


a P | 


-1/2 exp 




B 


\ 


Sl/=S m 

m 


3/2 


sinh 


/Xotii v T) 
V B ) 


(17) 


where u„, g„, and S*/ are, respectively, the reduced cyclotron mass, the effective 
g-factor, and the extremal cross-sectional area of the v sheet of the Fermi surface, 
Tq is the Dingle temperature representing the effect of electron scattering, and 
a = 27r^km o c/he. The phase factor is cp^ = 2n\y ±7t/ 4, where y is the phase associated 
with the Landau quantization; the upper sign is used when is a maximum, the lower 
for a minimum. From equation (15) we find the general expression for the strain: 

-Z4 

i=i 

6 

where the s?. are defined by > sP.c^ = 6., 
i] 9 iJ ik 

3=1 

enough so that the last term in equation (18) can be neglected. (Also, because of the 
smallness of M, the difference between B and H is inconsequential in our analysis.) 
The oscillatory magnetostriction in the [l , m, n] direction for a PbTe crystal with zero 
(or constant) external stresses is thus described by the equation 


— - 4ttM ^ 

a€ i 3€ i. 


(18) 




. In PbTe the magnetization is small 
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( 19 ) 


"Z mn 


\ ' 9 e i mn o 3f2 

L ■ ij ^ 


i,j=l 


where Be, /Be. is a product of direction cosines, 
L mn i r 


e Zmn = l2 H + m2e 2 + n % + mne 4 + lne 5 + Zme e 


The inverse of Young’s modulus in the [Z , m, n] direction is defined by 
-1 = 2 , /o„ ■ therefore the ME oscillations are given by 

6 




l mn' l mn 




"-E 


8€ Zmn s o 8g k 
9e i ^ 9cr Zmn 


i,j,k,t=l 


Be. 


x 


a 2 n 


aa k ae. ae t 


( 20 ) 


Using the identity 3or k / 9cr £ mn = 9e Znu/ 9e k and the a PP roximation 9e t/ 9cr k ~ s tk’ the 
ME oscillations are thus described to first order in the elastic constants by 


mn> 


6 

L 

i,j,k,t=l 


9e Zmn 9e Z mn O o O o a 2 fl 
S ik S jt — d€ t 


Be : 


8£ i 


ae>, 


( 21 ) 


Strain Derivatives of the Free Energy 

In order to evaluate the strain derivatives of the free energy S2 we need to evaluate 
the strain derivatives of the various band parameters which appear in fl. To this end, 
we next use our model of the PbTe band structure to express the strain dependence of the 
effective masses and cross-sectional areas in terms of the strain dependence of the 
Fermi energy and the band gap. We then obtain expressions for the strain dependence of 
these two characteristic energies in terms of the deformation potential parameters. 
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The band shifts described by the deformation potential scheme lead to variations in 
S2 via three types of effects. The first is due to the relative shift of the valleys pro- 
duced by a shear strain. Under a shear, some of the valleys move up in energy and 
some move down. This relative shift of the valleys is accompanied by intervalley charge 
transfer and changes in the Fermi surface cross-sectional areas. Second are effects due 
to the nonparabolic bands. Since the energy bands of PbTe are nonparabolic, a shift in 
the Fermi energy also produces a change in the band parameters. Third are the effects 
due to the strain dependence of the energy gap. The conduction and valance band defor- 
mation potentials are different; therefore, the energy gap is strain dependent. A change 
in the gap produces a change in the curvature of the bands, hence a change in band pa- 
rameters. The first effect can occur only with a shear strain, whereas the last two ef- 
fects will accompany both shear the dilatation strains. We will see from what follows 
that the dominant effect in PbTe is due to the intervalley charge transfer. This is evi- 
denced experimentally by the fact that both the [100] ME and MS oscillations, which are 
the only ones to come from pure dilatation strains and involve no intervalley charge 
transfer, are at least an order of magnitude smaller than the oscillations along other 
axes. 

We account quantitatively for the strain dependence of the cyclotron and band masses 
by differentiating equations (A9) and (A 5): 


where 


and 


illLM = c (m) + fi . c t ao] 3 - ln *** + C e (fi) — 

Be . 1 3e. L 1 J 3e. y 3e. 


3e. 

i 


3 In m. 


3 In m. 


= C ,(m*) I^L + [l - c t (m*)| 7 + C,(m*) 90 

ae i 1 d€ i L 1 J aej w 


3e i 


C t (tx) - 

w, 


cos 2 6 + sin 2 0 


2 V-i 


CJ/jl) = (-^} cos 6 sine (l - — 




H* 


(22) 


(23) 


(24) 


(25) 
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and C^(m*) and Cg(m*) are obtained from equations (24) and (25) by replacing p, p^., 
and p^ with m* , m^., and m^ , respectively. From equations (7) and (A9) we see that 
the transverse masses are functions of both the Fermi energy and the energy gap. There 
is a further energy -gap dependence of the band edge masses, m^ and m^, which we 
discussed in a previous section. Assuming that the effective gap is proportional to the 
direct gap, E <X A, simple differentiation gives the relations 

o 


3 In 


**t 


= C, 


d€. 


3 In E. 


9e i 
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3 In E 


g a* 


JL 


(26) 


and 


where 
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(27) 


(28) 


(29) 


and, further, 
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(30) 
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and 


where 


3 In m, 3 In E 
L g 


3e. 3e. 

i l 


C 

m 




(31) 


(32) 


Equations (22) to (32) express the strain dependence of p and m in terms of the 
strain dependence of the energy gap, the Fermi energy, and the angle between the 
spheroid axis and the magnetic field. 

From equation (A6) we see that the strain dependence of the extremal cross- 
sectional area is 


3 In S 3 In m* 3 In E„ 
m F 

S€j 3 £j atj 


(33) 


Thus, using equation (23), (30) and (30) with equation (33), we are also able to express 
3 In S m /3e. in terr.s of the strain dependence of E^, E F , and 6. 

We have assumed that E^, cc A; therefore 3 In E^/Se. = 3 In A/3e^, which is given 
by equation (13). The strain dependence of the Fermi energy is computed using the con- 
straint that the total number of carriers in the sample is independent of the strain. 
For our spheroidal model (see eq. (Al)), the density of carriers in the v in valley is 


P* 


,3/2 


3;r 2 h 3 


m t m Z 


1/2 



4 

Since £P = V ^ p^, where V is the crystal volume, the constraint 

v=l 

leads directly to 


(34) 


3^'/3e i 


= 0 
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= 0 (35) 


4 4 



v=l v=l 


where 6.. is the Kronecker symbol, and we have used equations (30) to (32) to express 
the strain derivatives of and in terms of the strain derivatives of and Ep. 
The Fermi energy for our p-type samples is the difference between the chemical poten- 
tial ? and the valance band maximum E^; therefore , 


9 In E^. _ t /9Ey 
8e i E F\ 9€ i 



(36) 


Using the valence band shifts and the gap dependence described by the deformation poten- 
tial scheme, equations (12), (13), (35), and (36) lead to expressions for 9 In in 

terms of the deformation potentials . There are two cases to be considered: 

(1) dilatational strain (i = 1, 2, or 3). Equations (12) and (13) indicate that all the 
gap and band-edge shifts are independent of v. Therefore, the sums in equation (35) 
are replaced by a multiplication factor of four, giving 


9 In E 


9e i 



(37) 


(2) shear strain (i = 4, 5, or 6). From equations (12) and (13) 


E 

v=l 
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From these and equations (35) and (36), it follows that 0£/0e. = 0, so that 


0 In E; 


3e. 





( 38 ) 


We next examine the strain derivatives of S2 which appear in the expression for the 
MS and ME oscillations, equations (19) and (21), using the expressions for the strain de- 
rivatives of the band parameters (eqs. (22) to (38)). Equations (16) and (17) indicate that 
dtl/de^ contains many strain derivatives since the effective masses, g-factor, Dingle 
temperature, cross-sectional area, and volume are, in principle, all strain dependent. 
However, by completely differentiating equations (16) and (17) and numerically evaluating 
the various terms using the parameters appropriate to our PbTe samples, we find that 
the contribution from the strain derivative of in the argument of the cosine factor is 
at least an order of magnitude larger than any other contribution. (In this numerical 

analysis, we assume that 9 1nT n /0e. and 9 In g/3e. are at most as large as 

^ ^ 1 _0 

0 In p/0^. Also, since the crystal strains never exceeded 10 in our experiments and 
the cubic symmetry was only slightly perturbed, we can assume that d6/de^ <1. This 
is equivalent to assuming that the Fermi surface ellipsoids continue to point very nearly 
along the body diagonals of the strained cubic lattice.) Schematically, 


= JA cos 

1 

8) 

1 

* , 

cS ™ 9 ln S m 

- A m m sin 

cS 

m - cp 


eKB 

efiB de . 

eKB 


Since we find that d In A/0e. is roughly of the same order as (or less than) 9 In s m / 9e j, 

one of the principal factors that determine the relative magnitude of the various terms 

is cS m /(efiB). This number is equal to 2ir times the number of Landau levels below the 

Fermi surface. Since our amplitude measurements were always made with oscillations 

corresponding to Landau numbers greater than five, the 0 In S m /0e^ term dominates. 

The fact that this is the dominant term is also borne out experimentally by the observed 

phase of the MS oscillations . The amplitude analysis in this work is made with that part 

of the data which show only a single period. These periods correspond to the smallest 

extremal area which we label S° . Since numerical evaluation shows that the contribu- 

m 

tions to 9S2/0e^ from other extremal areas and from harmonics (X > 1) are negligible, 
the expression that we use for our MS amplitude analysis takes the form 
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(39) 
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where the sum is over all the equivalent pieces of the Fermi surface which contribute to 
the given period. 

We further find from numerical evaluation of the various terms in the second strain 
derivatives of Q, that under the experimental conditions just discussed we can write 



Equation (40) could include terms in <^3^/0 e. 0e^ however, a comparison of phases of 
the MS and ME oscillations shows that such terms are negligible. 


Analytic Expressions and Ratio 

Equations (19) and (39) and equations (21) and (40) combine to give the expressions 
that we use in analyzing our MS and ME experimental data, 
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In our analysis of the experimental data, we fit the MS data to equation (41) and sepa- 
rately fit the ME data to equation (42). We determine experimentally the values of A 1Q 
and and find experimental values for the strain derivative sums which we then re- 
late to the deformation potentials. From equations (41) and (42) we see that the compli- 
cated coefficient A^q should disappear from the expression for the ratio of the MS and 
ME amplitudes for a given sample under the same thermodynamic conditions, that is, the 
same B, T, crystallographic orientation, etc. However, we find the MS and ME Dingle 

temperatures and to be different, leading to a term exp |ap J ej 

in the expression for the ratio. Even then, consideration of the ratio leads to a consid- 
erable simplification of the analysis, and we use the ratio in a separate approach to 
analyze the deformation potential. 


EXPERIMENTAL METHODS 

The oscillatory magnetostriction was measured by determining the change in capac- 
itance of a parallel -plate capacitor composed of one fixed plate and one plate glued di- 
rectly to the end of the oriented single crystal PbTe sample. We used a three -terminal 
capacitance cell based on a differential thermal expansion cell described by White 
(ref. 47), modified to allow a variety of samples with different dimensions to be used in 
one cell. A ratio transformer bridge (General Radio 161 5-A) and a phase -sensitive null 

detector were used to measure the change in capacitance. This technique allowed us to 

-fi 

detect a change in capacitance of the order of 10 picofarad, which corresponded to 
strains of a few times 10”^. 

The magnetoelastic oscillations were measured using the mechanical composite 
oscillator technique developed by Quimby and Balamuth (refs . 48 and 49) . The composite 
oscillator consisted of a quartz bar (typical dimensions 1 by 0.25 by 0.25 cm) bonded to 
an oriented single crystal PbTe bar of identical transverse dimensions. The long axis 
of the quartz transducer was in the Y quartz direction. The X-faces of the quartz were 
coated with chromium-gold electrode films so that a Young’s modulus vibrational mode 
for the composite system was excited when a sinusoidal voltage was applied to the elec- 
trodes. The electrical leads, which also served to support the oscillator, were attached 
to the transducer at nodal points. The composite oscillator was driven by a frequency- 
stable constant -voltage source, usually at a frequency near the fundamental resonance 
of the transducer (typically 250 kHz). Changes in the Young's modulus of the sample as 
a result of an applied magnetic field produced changes in the resonance frequency of the 
composite oscillator (appropriate analytic expressions are given by Balamuth (ref. 49)), 
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which, in turn, produced changes in the electrical impedance of the oscillator circuit. 

This impedance change was detected directly by monitoring- the voltage across a small 

series resistor. The voltage was amplified, rectified, and applied to an X-Y recorder. 

The response of the composite oscillator was investigated as a function of magnetic field 

at several driving frequencies on each side of the resonance as well as at the resonance 

peak itself in order to compare the changes in dispersion and attenuation of the system. 

In the 2 50 -kilohertz range, attenuation effects were completely negligible compared with 

the resonance frequency shifts due to changes in the PbTe Young's modulus. Because of 

3 4 

the high Q of the composite oscillator (typically 10 to 10 ), the sensitivity was high; 

— 

the minimum detectable fractional change in Young’s modulus was estimated to be 10" 

_7 

to 10 . It is noteworthy that even in its relatively unsophisticated form this system 

compares favorably in sensitivity with the "sing -around" ultrasonic technique commonly 
used to measure sound velocity changes (ref. 44). 

Details of the MS and ME experimental techniques and apparatus have been described 
elsewhere by one of us (T.E.T.) (ref. 50). 


EXPERIMENTAL RESULTS AND DISCUSSION 
Band Parameters 

Our experiments were performed on single -crystal samples from two different 
boules of p-type material grown by one of us (T.E.T.) by the Czochralski method. The 
boules were labeled X5 and X6. Samples were spark cut from neighboring sections of 
the middle of the boules. Hall measurements gave carrier concentrations of (3. 5±0.2) 
xlO per cubic centimeter for the X5 samples and (0.9±0.1)xl0 ° per cubic centimeter 
for the X6 samples. The X5[110] sample was a 0. 851- by 0. 488- by 0. 462- centimeter bar 
with the long axis parallel to [110]. The X5[100]-stack was 0. 927- by 0. 353- by 0. 36- 
centimeter [100]- oriented bar produced by gluing two nearly cubical pieces together with 
GE7031 varnish. The stacking was done to increase the OM signal. No magnetoelastic 
oscillation experiments were made with the samples from X5. The three X6 bars with 
axes along [100], [110], and [111] had dimensions 0. 749 by 0. 274 by 0. 312 centimeter, 

0. 745 by 0. 282 by 0. 244 centimeter, and 0. 805 by 0. 269 by 0. 257 centimeter, re- 
spectively. 

The oscillation periods and phases were analyzed in terms of the factors 
cos[ 27 t/(PB) - <pj]in equations (41) and (42). The period and phase were found by assign- 
ing successive integers N to the oscillation peaks and plotting the inverse magnetic field 
for the N** 1 peak, B^\ against N. The major oscillation periods found in this manner 
from the complete set of experiments are presented in table HI. The errors quoted are 
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TABLE III. - MEASURED OSCILLATION PERIODS AND CYCLOTRON MASSES 


Field direction 

Oscillation period, T * 

Cyclotron mass, m Q 



Sample 




X5 

X6 

X5 

X6 

Bll[100] 


(10. 9±0. 2)X10 -2 






BH[110] 

(5. 46±0. 03)xl0” 2 

(14. 6±0. 2) 

0.047*0.003 

0.033*0.003 

Bllflll] 

(6. 55*0.07) 

(18. 0±0. 3) 

.038*0.001 

.027*0.002 


a combination of the standard deviation of the mean found from a least-squares fit of the 
data and the 0. 5 percent field calibration uncertainty of the magnets. These periods are 
in agreement with those found in reference 27. 

An important consideration in the analysis was the relative phases of the MS and ME 
oscillations. As equations (40) and (41) show, the ME and MS oscillations for a particu- 
lar sample and field orientation should be 90° out of phase if our model is correct. We 
observe this phase relation in our data (fig. 3). The sign of the oscillations depends on 


o 

the sign of ( 3e Zmn/ 9€ P s ij( 9 ln S m //9e j^ and si £ n of -^io’ We tlia * a11 our 
i, j=l 

data fit consistently into our analytical picture if we assume that y = 1/2. This value 
of y leads to a negative value for Ajq, which from equation (17) corresponds to a nega- 
tive value for cos(?rgp/2). This sign for c os (7rgju/2) is in agreement with the gfj. val- 
ues presented in table I. A least-squares fit of all the phases within this scheme gives 
y = 0. 49±0. 01. 

For the spheroidal model described by equations (3), (5), and (6) the carrier den- 
sity p needed to fill four spheroids is related to the oscillation period P and the mass 
anisotropy K = m^/m^. by 


= aP" 3//2 (K cos 2 0 + 


sinV/V 1 / 4 


(43) 


where 8 is the angle between the major axis of the spheroid, that is, [111], and the 
applied magnetic field and a = 2.263X10 16 t" 3//2 cm" 3 (a product of universal con- 
stants). Using equation (43), the periods presented in table in, and the measured carrier 
concentrations, we find K = 6. 5(+l. 0, -0. 6) for X5 and K = 8. 3±0. 9 for X6. In figure 4 
we have plotted our values for K along with values found in p-type PbTe by other work- 
ers. The two most recent studies, those by Burke, Houston, and Savage (ref. 37) and by 
Schilz (ref. 32) indicate K values that are twice as large as those found in our work and 
that of three other investigators (refs. 26 and 51 to 53). Burke (ref. 37) made an extensive 
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Anisotropy, 



Figure 3. - Typical magnetostriction and magnetoelastic oscillation data. Sample X6 [110] ; field direction, 

B 1 1 [111]; Young's modulus and strain measured parallel to [1 10] direction; absolute temperature, 4. 2K. 
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Figure 4. - Published values of the (ill) cyclotron mass and the mass anisotropy K as function 
of carrier concentration for p-type lead telluride. 
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study of Shubnikov-de Haas frequencies. He Fourier -analyzed the frequency spectrum 
from the three fundamental oscillations and as many as five harmonics for numerous 
magnetic field directions in the (110) plane. It is quite difficult to understand the dis- 
crepancy between their results and the results of the other workers. It is noteworthy 
that Burke (ref. 37) finds that the carrier concentration determined from a high-field Hall 
measurement agrees with that determined from the best fit of their data to equation (43). 
In both the work of Stiles (ref. 52) and Cuff (ref. 26), the carrier concentration found 

from equation (43) was about 30 percent less than that determined from a Hall measure- 
18 

ment for p =* 3X10 per cubic centimeters. The latter authors attribute these extra 

carriers to other pieces of the Fermi surface. Analysis of Coste’s data (ref. 53) gives 

reasonable agreement between the p from equation (43) and that from a Hall measure- 
o 

ment. Since the value of K comes into the amplitude analysis through the term 
- 1/2 

\d^S u /dp^\ , (eqs. (17) and (A4)), the contradictory K values shown in figure 4 will 
be discussed further in the section Amplitudes. Unless otherwise stated we will use our 
experimental values for K. 

Cyclotron effective masses were found by measuring the ratio of the MS oscillation 
amplitudes at two temperatures. From equations (41) and (17) the ratio of the amplitudes 
at the same magnetic field but at two different temperatures is given by 


e(T l5 B) 
e(T 2 , B) 




* (44) 


This equation was solved by Newton’s method on a digital computer for a series of oscil- 
lation peaks and temperatures. Data were normally taken near 4.2, 6,5, and 8. 5 K. 
From a weighted least-squares fitting of the data result the masses presented in 
table HI. Since the MS amplitudes for Bil[100] were very small, no masses could be 
determined for this field direction. No temperature -dependence studies were made for 
the magnetoelastic oscillations. 

From equations (3) and (9), m^ = l/2(m^ + p 3 n i Q ) and m^ = Km^. Thus, from the 
values just presented, m t /m o = 0. 028±0. 002 and 0. 023±0. 002, and /m = 0. 18±0. 02 
and 0. 19±0. 02 for the X5 and X6 samples, respectively. 


%n his paper Coste claims that these two values do not agree, but his computation 
of p seems to be in error. 
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As a test of the internal consistency of the model, we computed the [110] cyclotron 
mass using equation (A9), m^, the measured [111] cyclotron mass, and the measured K. 
The computed values are 0. 044±0. 002 and 0. 032±0. 002 for X5 and X6, respectively. The 
agreement between these values and the measured values shown in table III supports the 
self-consistency of the model and measurements. 

The Fermi energy is easily computed from the de Haas-van Alphen period and the 
effective mass (see the appendix): 


E 


F “ 



(45) 


From our data for Bllflll], where m* = m^., we find Ep = 63±5 and 28±2 meV for 
samples X5 and X6, respectively. 

Above 2 tesla (20 kG) the ME and MS oscillations begin to show a splitting of the 
oscillation peaks due to the Zeeman splitting of the Landau levels. From the oscilla- 
tions for Bll[100], where all four ellipsoids have the same fundamental frequency, we 
were able to determine the spin splitting of the peaks. Each peak corresponds to a solu- 
tion of the equation 


J^ = N±M +y+ l (46) 

PB 4 8 

where the second term on the right accounts for the two spin levels for each Landau 
quantum number N. This term was found by indexing the oscillation peaks with a spin 
index and a quantum number and plotting N against (B^) + and (B^)_. The resulting 
g-factors are compared in table I with those found in p-type material by other workers. 
Butler and Calawa (ref. 34) obtain their results from magnetoemission sutdies in diode 
lasers, Schilz (ref. 32) from magnetoacoustic attenuation oscillations, and Cuff (ref. 27) 
and Burke (ref. 37) from Shubnikov-de Haas oscillations. For both Bllflll] and Bll[110], 
there are additional oscillations from a second set of ellipsoids. This makes the analysis 
difficult, and it has not been carried through. 

The value of gp depends on the assignment of spin and quantum numbers to a set of 
oscillations. From the phase analysis mentioned previously, we know that cos(xrgp/2) 
must be negative, but this condition does not determine a unique labeling of the oscillation 
peaks. The values for gp/2 in table I correspond to a spin splitting greater than half 
the Landau splitting. This scheme gives g values that agree with the values found by 
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Patel and Slusher (ref. 54) from spin-flip Raman scattering in n-type material with 
17 

n = 10 per cubic centimeter. 


Amplitudes 


The MS amplitude coefficients were found by fitting the strain oscillation amplitudes 
to equations (17) and (41). From the data we computed the quantities y = lnje^ mn B -1//2 

x sinh(o'pT/B)| for each oscillation peak and a linear least-squares fit was made to a 

1 

plot of y against B . From the equations, it follows that the slope of this plot gives 
the Dingle temperature and that the y intercept leads to a value for 
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This procedure necessitates a knowledge of S^, p, and 


) 2 S/Sp£ 


- 1/2 


S=S L 


for which we 


m 


use the band parameters presented in the section Band Parameters and equation (A3). 
The Dingle temperatures found from the MS oscillations are presented in table IV. The 


TABLE IV. - DINGLE TEMPERATURES FROM 
MAGNETOSTRICTION OSCILLATIONS 


Sample 

Strain 

Field direction 


direction 

Bll[110] 

BH[110] 

Bll[lll] 

Bll[lll] 



Dingle temperature, 

K 

X5 

[110] 

5. 4±0. 5 

5.2±0. 5 


4. 5±0. 3 


X6 

[110] 

9 . 2± 0 . 4 

8. 3±0. 5 

7. 2±0. 3 

8. 3±0. 6 


[111] 


7. 3±0. 5 






^Because of the similarities between the PbTe conduction and valance bands, the 
g-values (as well as effective masses) should be comparable in n- and p-type material. 
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errors indicated are the standard deviation of the mean computed from the random er- 
ror associated with the amplitude measurements and the uncertainty in p. The mag- 
netostriction amplitudes for X5[100] and X6[100] were too small to give reliable Dingle 
temperatures. 


In table V the strain derivative sums determined from the experimental MS ampli- 
tudes are presented and compared with values computed from the theory presented in the 
section Strain Derivatives of the Free Energy. The theoretical values were computed using 
the deformation potentials from Ferreira’s (ref. 41) theoretical work (table II), and the 
elastic constants measured by Houston, Strakna, and Belson. 4 The extrapolated 
T = 0 K values for the stiffness constants were used to compute the following compliance 
constants: s° 4 + 2s° 2 = (0. 731±0. 007)xl0 -11 m 2 N -1 and s^ 4 = (6. 61±0. 04)xl0 -11 
m^N -4 . The large uncertainties of the theoretical values result entirely from the un- 
certainty in cos(jrgp/2). Even though gp is experimentally determined to better than 
10 percent, cos(ngp /2) has a large percentage uncertainty when the function nears zero. 
The values in table I lead to cos(ffgp/2) = -(0. 19±0. 10) for X5 and -(0. 40±0.08) for X6. 
The errors in the experimental strain derivative sums are the standard deviation of the 
mean due to the uncertainty of the amplitude measurements and the uncertainties in p 

and K. Because of the discrepancy between the K values found by different experi- 

18 

menters for p s=s 3x10 per cubic centimeter there is an added degree of uncertainty 
about the experimental values for X5, since the amplitude analysis includes the term 


a 2 sV3p 


1 - 1/2 


which is K dependent. For instance, if we were to use K = 13, the 
value of Burke, Houston, and Savage (ref. 40), rather than our experimental value 
(K = 6.5) all of the X5 experimental values in table V would be reduced by 23 percent for 
Bll[100], by 27 percent for Bll[110], and by 29 percent for B ll[l 1 1 ]. Using K = 13 in- 
stead of 8.3 for X6 gives 17-, 19-, and 20-percent reductions, respectively, for these 

o v 2 - 1 / 2 

three cases. Since both cos(jrgp/2) and 9 S /9pg cancel out of the expressions 

for the ratio of the ME and MS amplitudes, the large uncertainties of these parameters 
do not enter into the results of the ratio analysis. 

Magnetostriction oscillations were investigated at temperatures from 4 to 9 K. No 
temperature dependence of the strain derivative sum was observed and the data in table V 
are the average of values over this temperature range. The oscillations in several 
cases were just barely detectable or were not seen at all. For these cases an upper 
limit for the amplitude coefficient is presented. We had no X5[lll] sample. Because 
the amplitudes for X6[100] were quite small for transverse magnetic field configurations, 


Straight forward, although tedius, symmetry considerations appropriate to the 
m3m symmetry of PbTe can also be used to find these relations. 
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TABLE V. - EXPERIMENTAL AND THEORETICAL VALUES OF STRAIN DERIVATIVE SUMS a FROM MAGNETOSTRICTION OSCILLATIONS FOR SAMPLE X6 (p = 9xl0 17 cm" 3 ) 

AND SAMPLE X5 (p = 3. 5xl0 18 cm" 3 ) 

[Theoretical values were computed using 2^ = 10. 46 eV. ] 


Sample 

Type of data 

Strain direction 

b [100] 

(110] 

fin] 

Field direction 

[1001 

[010] 

[011] 

[110] 

[110] 

[001] 

[111] 

[111] 

[110] 


Strain derivative sum a , m 3 N” X 

X5 Experimental <0.97x10" 10 

Simple theory (0. 04*0. 02)xl0' 10 
Full theory (. 19*0. 10) 

<0.42x10" 10 
(0. 04±0. 02)xl0" 10 
(. 19±0. 10) 

<0.34x10" 10 
(0.02±0.01)xl0“ 10 
(. 13*0. 07) 

(6.5*0. 5)xl0" 10 
(6. 9*3. 6)xl0" 10 
(9. 2*4. 9) 

(6.3±0.7)xl0“ 10 
(7.0±3.7)xl0" l ° 
(9. 5*5.0) 

<0.85xl0" 10 
(0. 04*0. 02)xl0~ 10 
(. 19*0. 10) 

(3. 5*0. 4)X10" 10 
(3. 5*1. 8)xl0” 10 
(4. 8*2. 5) 

(c) 

(6.9±3.6)xl0~ 10 
(9. 4*4. 9) 

(O 

(4. 6*2. 4)X10~ 10 
(6. 4*3. 4) 

X6 Experimental <4.7 v 10~ 10 

Simple theory (0. 08*0. 02}xl0~ XX) 
. Full theory (.30*0.06) 

(0.08±0.02)xl0" l ° 
(. 30±0. 06) 

{0.04*0. 01)xl0" 10 
(. 18*0.04) 

(50±4)xl0" 10 

(32.9*6.6) 

(38.9*7.9) 

(56±9)xl0~ 10 

(32.9*6.6) 

(39.3*7.9) 

<5. 7x10" 10 

(0. 08*0. 02)xl0" 10 
(.30*0.06) 

(22±l)xl0“ 10 

(16.5*3.3) 

(19.8*4.0) 

(38±3)xl0‘ 10 

(32.9*6.6) 

(39.2*7.8) 

(22±3)X10’ 10 

(22.0*4.4) 

(26.2*5.2) 

a Strain derivative sum = 

c-g,'2£ J ln S m ?f j> 

s° 



b Sample X5 was stacked to increase the MS signal. 
c No sample. 



we did not attempt to aline the field along specific crystallographic directions for trans- 
verse fields. These facts explain the absence of entries for the corresponding crystal 
and field orientations in table V. 

Two sets of theoretical values are presented in table V. The "simple theory" 
values correspond to a rigid parabolic band model with infinite energy gap. In this model 
the strain mass term, 3 In m* /3e^ in equation (33), is neglected and 3 In S^/Se^ is 
considered to be due solely to the change in crystal volume (= -2/3) for dilatation strains 
(i = 1, 2, or 3) or due to intervalley charge transfer (= 1/3 E^/E^j for shear strains 
(i = 4, 5, or 6). The values labeled "full theory” are computed using the model repre- 
sented by equations (33), (22), (30), (31), (37), (38), and (13). As an example of the 
magnitudes of some of the numbers which enter the theoretical computations, 3 In m* 
/3e., 3 In E-p/3^, and 3 In S^/Se^ for X6 with Bll[110] are presented in table VI. 


TABLE VI. - REPRESENTATIVE THEORETICAL STRAIN 
DERIVATIVES FOR X6 WITH Bli[100] 


Type of theory 

Value of i 

Strain derivatives 

3 In m* 

3e i 

3 In Ep 

ainS m 

Simple 

1, 2, or 3 

0 

-2/3 

-2/3 


4, 5, or 6 

0 

124.5 

124. 5 

Full 

1, 2, or 3 

12.6 

-15.7 

-3. 1 


4, 5, or 6 

23. 5 

124. 5 

148.0 


The agreement between the experimental and theoretical values in table V is quite 
reasonable; however, due to the large uncertainty in the theoretical values, especially 
for the X5 data, a definitive quantitative analysis of the experimental and theoretical am- 
plitude factors is not possible with the numbers presented in table V. The relative mag- 
nitudes and ratios of the various terms do confirm many of the aspects of the model we 
have presented. We can see this by writing the strain derivatives in a more transparent 
form. When evaluating the Strain derivative sums for the strain and field directions 
shown in table V by means of the prescriptions described in the section Strain Deriva- 
tives of the Free Energy, one finds that, except for the case el1[110] with Bll[001], all 
the sums can be simplified into the sum of a dilatation term and a shear term of the form 
(ref. 55) 
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9 In S’* 

s o m = Ne( ^ 


i] 


9e 


+ D/3 g ) 


(47) 


where N g is the number of equivalent extremal areas S m’ % " ( s ll + 2s 12J d In 

/3e, , 0 = s?. 0 In S° / 0 €., and D is a direction cosine factor which takes the values 0, 
1/2, -1/2, 1, and -1/3 for the directions (100), [110], [110], [111], and [111], respec- 
tively, when the derivative indicated in /3 g is performed on the extremal area of the 
[lll]-spheroid. For the case ell[110] with Bll[001], the right side of equation (47) be- 
comes simply 4/3^. A significant aspect of the model is that in this case the shear terms 
cancel when summed over the four spheroids: the shear strain causes two of the valance 
band maxima to move up in energy, thereby increasing the extremal areas at the Fermi 
energy, and the other two move down, decreasing their areas. Comparing equation (47) 
with the experimental data in table V shows that, qualitatively, the data fit our theoretical 

model quite well. We find, for example, that I 6 I » I B, I; the full theory values in 

-11 s a -11 
table VI lead to j3^ = -2. 2x10 square meters per newton and /3 = 977x10 square 

meters per newton. These values are typical of the large difference between the dilata- 
tion and shear contributions to the strain derivative sums predicted by our theoretical 
model, including other field directions and for sample X5. Since the shear term /3 g is 
theoretically so much larger than the dilatation term /3 d , any term in equation (47) that 

contains R will be much larger than terms which do not contain fl. 

s s 

From table V we see that for all cases where equation (47) shows contributions from 
only the dilatation strains (i.e. , ell[100], D = 0, or ell[110] with Bll[001], the experi- 
mental values are at least an order of magnitude smaller than those for which equa- 
tion (47) shows shear contributions. This result confirms one aspect of the model. 

Another confirming aspect of the MS data comes from the ratios of various ampli- 
tude terms. Using equation (47) and assuming_that /3 g » (3^, the strain derivative sums 
for gll[110] with e along [110], [110], _and [111] should be in the ratio l/2:-l/2:-l/3. 

For Bil[lll]with e along [111] and [110] the ratio should be 1: -1/2. The data, as pre- 
sented in table V and interpreted in light of the symmetries involved, are seen to be in 
good agreement with these predictions. Perhaps the most striking experimental con- 
firmation of the model comes from the signs of the observed phases. As we pointed out 
the phases are in complete agreement with our model; for example, the MS oscillations 
for BH[110] with ell[II0] and [110] are 180° out of phase, as we would expect from equa- 
tion (47). 
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The observed changes in amplitude and the phase shifts of the MS oscillations due to 
cancellation of the shear terms as the magnetic field is rotated into the [001] direction 
in the (110) plane are shown by the experimental data in figure 5. The experimental re- 
sults with B or e along (100) are difficult to analyze because the strain amplitudes 
are quite small and any slight misalinement of the sample would produce an anomalously 
large signal. This fact is made quite apparent by figure 5. 


180 


120 


60 


0 <— 



Figure 5. - Amplitude and phase of [no] magnetostriction oscillations for B near 
[003 for sample X6. 


The ME amplitude analysis is quite similar to the MS analysis. A linear least 

squares fit was made of x = ln|A(Y^ mn ) -1 B 1//2 sinh(o!fji.T/B)J against B" 1 . From 

equations (17) and (42) it follows that the Dingle temperature comes from the slope of the 
curve and that the x intercept gives a value for 



The Dingle temperatures found from the ME amplitude analysis are presented in 
table VII. The ME Dingle temperatures for Bll(llO) are consistently lower than the MS 
Dingle temperatures for the same samples, whereas for Bll(lll) the ME Dingle tem- 
peratures are larger than the MS Dingle temperatures. This result is unexpected within 
the framework of our model. 
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TABLE VH. - DINGLE TEMPERATURE FROM MAGNETOELASTIC 


OSCILLATIONS 


Sample 

Field direction 

and axis 

Bll[100] 

Bll[001] 

BN[110] 

B!l[110] 

Bllflll] 

Bllflll] 



Dingle temperature, 

K 


X6[100] 

9. 6±0. 8 

8 . 3± 1 . 0 


6 . 8± 0 . 9 






X6[110] 

X6[lll] 


8. 4±0. 4 

8. 5±0. 4 

7. 5±0. 3 
7. 3±0. 5 

10. 2±0. 2 

8. 2±0. 5 






Table Vm shows theoretical and experimental values determined for the ME strain 
derivative sums. The simple theory and full theory values have the same meaning as 
the MS theoretical values. The elastic constants, deformation potentials and the value 
of cos(7rgp/2) used for the theoretical evaluations are the same as those used to calcu- 
late the MS amplitude terms. Again, the uncertainty in the theoretical values comes 
from the uncertainty in cos(7rgp/2). Recall also that the experimental values depend on 
the value of K used in the data reduction. Using K = 13 instead of 8. 3 leads to a re- 
duction of the experimental values by 17, 19, and 20 percent for the [100], [110], and 
[111] magnetic field directions, respectively. 

Just as we were able to write the MS strain derivative sums in a form that illumi- 
nates the signficance of the relative magnitudes of the values in table V, we find that the 
ME strain derivatuve sums for all the sample and field directions shown in table VII ex- 
cept the case e M[1 10] with Bll[100], can be expressed in the form 



where all the terms have the same meaning as in equation (47). For the case ell[110] 
with BilflOO], the right side of equation (47) becomes 4^ + ^ /3 g ). Remember that, 
for the comparable MS case, the shear terms R cancelled out when summed over the 
four spheroids. In the ME case, the strain derivatives are all squared before they are 
summed, and the shear term does not cancel. Since |/3 g | » 1/3^1, as seen experi- 
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TABLE VIII. - EXPERIMENTAL AND THEORETICAL VALUES OF STRAIN DERIVATIVE SUM FROM MAGNETOELASTIC 


OSCILLATIONS FOR SAMPLE X6 
[Theoretical values were computed using = 10. 46 eV. ] 


Type of data 


Sample and axis 


X6[100] 


X6[ 1 10 ] 


X6[ 111] 


Field direction 



[100] 

[010] 

[Oil] 

[110] 

[110] 

[001] 

[111] 

[111] 

[110] 

1 

1 




Strain derivative sum a , 

m 4 N' 2 




Experimental 

(3.1±0.7)xt0“ 18 

(2.S±0.9)xl0" 18 

(1.3±0.3)xl0“ 18 

(23±3)xl0‘ 18 

(15±2)xl0' 18 

(33±4)xl0‘ 18 

(9.7±1.6)xl0‘ 18 

(2.7±0.4)xl0’ 18 

(2.6±0.3)xl0* 

Simple theory 

0 

0 

0 

(13.5±3) 

(13.5±3) 

(27±5) 

(6.8±1.4) 

(27±5) 

(6±1) 

Full theory 

0 

0 

0 

(18.9±4) 

(19 ,3±4) 

(37±7) 

(9.8±2.0) 

(38.5±8) 

(8.6±1.7) 


a Strain derivative sum = | cos(irgp/2) 


v 


6 

E 






mentally and theoretically from both the MS and ME data, the MS strain derivative sum 
for the case eli[110] with Bll[100] was an order of magnitude less than the sums in- 
volving ft . In contrast, the ME strain derivative sum for this orientation is of the same 

b 

magnitude as in the other cases. 

Using |jS g | » 1/3^1, equation (48) predicts that for X6[110] the ratio of the strain 
derivative sums for Bll[110], [110], [001], and [111], respectively, should be 2: 2: 4:1 
(table VHI). The experimental values confirm this aspect of the theoretical model rea- 
sonably well. The experimental values for X6[110] are also in reasonably good agree- 
ment with those predicted by the theoretical model. The fact that the ME theoretical 
values predicted for X6[100] are much smaller than the values measured is not too sur- 
prising since our model and calculations contain several simplifying assumptions. The 
X6[100] experimental values are, in fact, an order of magnitude smaller than the 
X6[110] experimental values. This we understand qualitatively within the framework of 
our model. The order of magnitude disagreement between experiment and theory for 
X6[lll] is not understood at this time. The experiments on X6[lll] were carefully re- 
peated in an unsuccessful attempt to resolve this discrepancy. Equation (48) predicts 
that the strain derivative sum for X6[lll] with B ll[ 1 1 1 ] should be equal to the X6[110] 
value for Bii[001] and that the value for the X6[lll] case with Bll[110] should be 2/9 of 
this value. The MS data for X6[lll] gave good agreement between experiment and 
theory, so the discrepancy does not seem to be associated with this particular sample. 

We have already pointed out that the expression for the ratio of the ME to the MS 
amplitudes is greatly simplified because of the cancellation of most of the coefficient 
Ajq which contains, among other things as factors, K and cos( 77 gp/ 2 ). From equa- 
tions (41), (42), (47), and (48), we see that, for all the sample configurations in 
tables V and VII except the case ell[110] with Bll[001], the expression for the ratio of 
the amplitudes takes the form 


A(Y 7 ) 

L mn 


-1 


imn 




d€ 


s ° il nS ™ exp 
3f i 11 8£ j 


a/i 



(49) 


If the Dingle temperatures were the same the right side of equation (49) would have a 
simple B~* field dependence. The experimentally observed field dependence does, in 
fact, deviate from B”* and we attribute this wholly to the difference in Dingle tempera- 
tures. 
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Because the X6[100] amplitudes are too small and the X6[lll] ME data are ques- 
tionable, we have only the data from X6[110] for the ratio. (No ME measurements were 
made on the X5 samples.) By least-squares fitting a semilogarithmic plot of equa- 
tion (49) for the data from X6[110], we find the parameters presented in tablex IX and X. 
As can readily be seen from table IX, the ME and MS Dingle temperature difference de- 
termined from the ratio (eq. (49)) and the difference obtained from the values from the 
two separate experiments agree. The nature of this difference is not understood within 
the framework of our thermodynamic model. 

The strain derivative sums determined from the ratio are displayed in table X. The 
terms simple theory and full theory have the same meanings as above. The experi- 
mental values for the ratio are independent of the K value and the cos(jrgju./2) so that 
uncertainties in these parameters have no effect on the results. The uncertainties of the 
experimental values in table X are mainly due to the random uncertainties of the ampli- 
tude measurements themselves. 

From our theoretical model (see, for instance, table VI and eq. (38)) we find that 
the strain derivative sums shown in table X are essentially linearly proportional to E^; 
therefore, for the theoretical values in table X we have scaled the deformation potential 
to give the best fit of the theory to the experiment. The best value EX = 7- 9±1. 3 elec- 
tron volts closely agrees with the 8. 5- electron-volt value found from piezo resis- 
tance measurements by Burke (ref. 51). When we scale the ME data for a best fit be- 
tween experiment and theory, we find E^ = 10. 1±1.2 electron volts. The uncertainty in 
the deformation potential determined from the ME data comes mainly from the un- 
certainty in cos(3rg|u/2). For the ME case, however, we must concern ourselves with 
the question of whether K is 13 or 8. 3. The 10-percent uncertainty introduced by this 
difference is not included in the ±1.2 electron volts. 


TABLE IX. - THE DIFFERENCE BETWEEN MAGNETOSTRICTION 
AND MAGNETOELASTIC DINGLE TEMPERATURES 
FOR X6( 110] 


Values from - 

Field direction 


[110] 

(U0| 

fiiil 


€ Y 

Dingle temperature difference, T^ - T^, K 

Eq. (49) a 

l. 4±0. 1 

1. 3±0. 7 

CM 

o 

1 

Tables IV and VII b 

. 7±0. 6 

. 8±0. 6 

- 1 . 0± 0 . 6 


a Fitting the ME-MS ratio data to eq. (49). 

b The difference between the values determined by the two 


separate experiments. 
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TABLE X. - STRAIN DERIVATIVE SUM DETERMINED 


FROM MAGNETOSTRICTION -MAGNETOELASTIC 
AMPLITUDE RATIO FOR X6[110] 

[Theoretical values computed using E v = 10. 46 eV. ] 


Type of data 

Field direction 


[110] 

[110] 

[111] 


Strain derivative sum 3 

, m^' 1 

Experimental 

(37±2)X10~ 10 

(26±5)xlO' 10 

(48±3)xl0" 10 

Simple theory 

(31±5)xl0~ 10 

(31±5)xl0' 10 

(31±5)X10' 10 

Full theory 

(37±6)xl0~ 10 

(37±6)xl0“ 10 

(37±6)X10' 10 


Strain derivative sum = 


L < S W 3 c Xj< 31 nS m/ a€ j> 

i,3=l 


TABLE XI. - SPIN SPLITTING FACTOR (ABSOLUTE VALUE) DETERMINED 
FROM MAGNETOSTRICTION AND MAGNETOELASTIC AMPLITUDE DATA 


]H^ - 7. 9*1.3 eV. 



Sample 

Field direction j 


and strain 
direction 

[110] 

[110] 

[0011 

[1111 

[111] 




Spin splitting factor, | cos(7igp/2) | 


Magneto- 

striction 

X 5[ 1 10 ] 
X6[110] 
X6[ 111 1 

0. 18*0.03 
.68*0. 13 

0. 17*0.03 
.76*0. 18 
.44*0. 10 


0. 18*0.04 
. 59*0. 10 



0. 52*0.10 





Magneto - 
elastic 

X6[ 110J 

0.86*0.33 

0. 55*0. 21 

0.63*0.24 

0.70*0.28 




Having found a value for the deformation potential from the MS -ME amplitude ratio 
we next return to the MS and ME data to find a value for cos(7rgp/2). Taking = 

7. 9±1. 3 electron volts, we use our theoretical model to compute values for the strain 
derivative sums and then combine these results with equations (41) and (42) and the ex- 
perimental amplitudes to find cos(77gp/2). The results of this procedure are given in 
table XI. The comparison of the spin factor found in this manner with the value found 
directly from the splitting of the oscillation peaks is a further test of the internal con- 
sistency of the theoretical model. Assuming cos(jrgp/2) to be the same in the [100], 
[110], and [111] directions, and incorporating our phase results, which indicate that this 
function must be negative, the average of the X5 values in table XI give cos(7rgp /2) = 

-(0. 18±0. 02) compared with the spinsplitting value of -(0. 19±0. 10). The X6 values 
average to -(0. 57±0. 06) compared with the spin splitting value of -(0. 40±0. 08). For this 
analysis we have used our own values for K. If we had used the K value of Burke, 
Houston, and Savage (ref. 40) the X5 and X6 data in table XI would have to be shifted 
down by 28 and 15 percent, respectively. If our model is correct, this indicates that a 
value of K = 13 for p = 3x10* is too high. 


CONCLUSIONS 

A detailed experimental and theoretical study of quantum oscillations in the mag- 
netostriction (oscillatory magnetostriction) and Young's modulus (magnetoelastic oscil- 
lations) of p-type lead telluride has been presented. The Young's modulus oscillations 
represents the first observation of quantum oscillations in the elastic constants of a 
semiconductor. In the theory, the oscillatory magnetostriction and magnetoelastic oscil- 
lations are related to strain derivatives of extremal Fermi surface cross-sectional areas 
using appropriate thermodynamic derivatives of the modified Lifshitz-Kosevich expres- 
sion for the quantum oscillatory electronic free energy. These strain derivatives are, 
in turn, related to a deformation potential description of the strain dependence of the 
electronic band structure using an adaptation of Cohen's nonellipsoidal , nonparabolic 
band model to describe the nearly spheroidal nonparabolic valence band of lead telluride. 

From the periods, phases, and spin splitting of the oscillations, we find experi- 
mental band parameters for p-type lead telluride that are in generally good agreement 
with parameters reported by other workers, although there is a discrepancy with recent 
results of Burke, Houston, and Savage (ref. 32) and of Schilz (ref. 37) in the value of the 

Fermi surface anisotropy constant K. These workers found values of K near 13 for 

18 

carrier concentrations near 3x10 per cubic centimeter, whereas we found K = 6. 5 
18 

for p = 3. 5x10 per cubic centimeter. 
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The absolute amplitudes of the magnetostriction and magnetoelastic oscillations were 
were generally in good agreement with the predictions of our theory. From the ratio of 
the amplitudes of the two types of oscillation, we have derived a value of the valance- 
band deformation potential = 7.9±1. 3 electron volts. This agrees well with a value 
previously obtained by Burke (ref. 16) from piezoresistance measurements. The rela- 
tive phase of the two types of oscillation is in accord with that predicted by the theory. 

A detailed study of the absolute amplitudes of the two types of oscillation indicates that 
the dominant effect of strain on the valance band is a lifting of the degeneracy of the 
valence band maxima by shear strains accompanied by intervalley charge transfer. The 
effects of strain on the band parameters (the so-called strain-mass effects) are roughtly 
an order of magnitude smaller than the intervalley charge transfer effects. This is in 
contradiction with a conclusion of reference 16 based on oscillatory magnetostriction 
experiments only. In that analysis these workers used band parameters from the litera- 
ture that give a value of the spin-splitting factor in disagreement with our experimental 
observations, and, in addition, their analysis appears to contain a numerical error. If 
these two features of their analysis are corrected, their results are consistent with 
small strain-mass effects. Study of our amplitude results also suggests that our value 
of the Fermi surface anisotropy constant is more nearly correct than the higher values 
reported by Burke, Houston, and Savage (ref. 37) and by Schilz (ref. 32). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 6, 1971, 

112-27. 


38 



APPENDIX - SOME BASIC RELATIONS 


In this appendix we present some basic relations between the band parameters for 
our spheroidal, nonparabolic energy band model of PbTe. Most of the relations follow 
directly from geometric considerations. 

The four (111) spheroidal bands described by equations (3), (5), and (7) can accom- 
modate a total carrier density p if the volume V p in momentum space of each sphe- 
roid is 


V e f> 


/E\h 3 = 1 ir2m t E F (2m ; E F ) 1/ ' 2 


This requires that the Fermi energy Eg be 


IT 


37T p 
4 


2/3 


E F = 


2m^ K 


1/3 


(Al) 


From straightforward geometrical considerations, the area of intersection S be- 
tween a plane perpendicular to an axis Pg and a spheroidal Fermi surface is found to be 


S ( 0, Pg) — 


27rm t Vk Ej 




2 V^ 


(k cos 2 0 + sin 2 0^ cos 2 8 + sin 2 0^ 


3/2 


(A 2) 


where 9 is the angle between the Pg axis and the spheroid axis and Pg is the dis- 
tance from the plane to the center of the spheroid. 

The period of the quantum oscillations is determined by the maximum cross- 
sectional area S m of che Fermi surface perpendicular to the field direction, 
p = eh/(cS m ). From equation (A2) we see that the maximum cross-sectional area occurs 
when Pg = 0, giving 


He) = 


eK^K cos 2 0 + sin 2 0^ 

H^ e f 


cm. 


(A3) 


From equation (A2) we can evaluate one of the important amplitude factors of the 
Lifshitz-Kosevich equation (eq. (17)): 
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a 2 s 


9p | 


- 1/2 


k / /2 Lsh ♦ =to 2 e \ S/4 


27T/ 


K 


(A 4) 


S=S. 


m 


It is convenient to define a band mass m* for the nonparabolic bands by the rela- 


tion 


J_\ 2 s cos 2 e + sin 2 0 


(A 5) 


m. 


l t 


The preceding expressions can be written more easily in terms of this band mass: 


S m (6) = 27rai*E F 


(A 6) 


P(0) 


eh 


cEp,m 


(A 7) 


a 2 s 


0p B 


- 1/2 


27 r/ \m*/ 


3/2 


(AS) 


S=S^ 


m 


The reduced cyclotron mass for a magnetic field in the B-direction is defined as 

ju = (l/2irm o ) ( 9S m / 9E )E • Substituting equation (7) into equation (A2) and differentiating 

F 

result in 


where 



co s 2 6> + sin 2 0 1 
l 2 HV-i 


$ 


(A9) 
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and 


m 3 



$ = (1 - £)(1 + 2r) 2 (l + 2r - rw)' 2 


and where 


and 



on = - (1 + r) 
2 


+ r + 


^ co t 2 ej 

m i / 


ij = r (l + 2r + — - cot 2 0 

\ m, ) 


Using the experimental results for PbTe <j> is found to be very near 1 when the field is 
in any of the principal c ry stal 1 og r aphi c directions: For the X6 samples 0.99 < $ < 1, 
and, for X5, 0.97 < < 1. Thus, for all practical purposes we can set $ = 1. Notice 

that / m^/m Q , except for Ep/E^ = 0. 
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